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Abstract 

We study the temperature dependence of the monopole condensate in different Abelian 
projections of the SU{2) lattice gauge theory. Using the Frohlich-Marchetti monopole 
creation operator we show numerically that the monopole condensate depends on the 
choice of the Abelian projection. Contrary to the claims in the current literature we 
observe that in the Abelian Polyakov gauge and in the field strength gauge the monopole 
condensate does not vanish at the critical temperature and thus is not an order parameter. 
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,1 Introduction 

m . 

O -The confinement of color in QCD is one of the most interesting issues in the modern quantum 
field theory. Numerical simulations of non-Abelian gauge theories on the lattice pQ show that 
the confinement of quarks happens due to a formation of the chromoelectric strings spanned 
| between quarks and anti-quarks. Although an analytical derivation of the color confinement is 
not available yet, there are various effective models which describe the emergence of the string. 
According to the dual superconductor model [2J, the vacuum of a non-Abelian gauge model may 
be regarded as a medium of condensed Abelian monopoles. The monopole condensate squeezes 
the chromoelectric flux (coming from the quarks) into a flux tube due to the dual Meissner 
effect. This flux tube is an analogue of the Abrikosov vortex in an ordinary superconductor. 

The basic element of the dual superconductor picture is the Abelian monopole. This ob- 
ject does not exist in QCD on the classical level, but is can be identified with a particular 
configuration of the gluon fields with the help of the so-called Abelian projection The 
Abelian projection uses a partial gauge fixing of the SU(N) gauge symmetry up to an Abelian 
subgroup. The compactness of the residual Abelian subgroup guarantees the existence of the 
Abelian monopoles in the Abelian projection. 

Many numerical simulations show that the Abelian degrees of freedom in an Abelian projec- 
tion are responsible for the confinement of quarks (for a review, see, e.g., Refs. j3]). One of the 
striking features of the Abelian projection is the effect of the Abelian dominance the Abelian 
gauge fields provide a dominant contribution to the tension of the confining string. Moreover, 
the internal structure of the string energy, such as energy profile and the field distribution are 
very well described by the dual superconductor model |T]. 

Since qualitative features of the confinement mechanism in QCD and in the pure £77(2) 
gauge theory are expected to be the same, below we restrict ourselves to the simplest case 
of the £77(2) gluodynamics. Most of the results supporting the dual superconductor scenario 
were obtained in the so called Maximal Abelian (MA) projection According to numerical 
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simulations [S] the monopole condensate in the MA gauge is formed in the low temperature 
(confinement) phase and the condensate disappears in the high temperature (deconfinement) 
phase in the perfect agreement with the expectations coming from the dual superconductor 
scenario. 

Besides the MA projection there are Abelian projections which are defined by a diagonal- 
ization of certain adjoint operators X[U] with respect to the SU(2) gauge transformations [3]. 
After the Abelian projection is fixed, the matrix X[U] becomes diagonal and the theory pos- 
sesses the (residual) U(l) gauge symmetry. The most popular examples of such gauges are the 
Abelian Polyakov (AP) gauge and the Abelian field strength gauge (F 12 gauge). These gauges 
correspond, respectively, to the diagonalization of the Polyakov line and the Uu plaquette 
variable. 

One may ask whether a dual superconductor nature of non-Abelian vacuum is realized in 
all Abelian gauges. Needless to say that this question is important for understanding of the 
properties of the QCD vacuum. Indeed, it seems natural that the color confinement must 
be described by a projection-independent model since the confinement is a gauge invariant 
phenomenon. On the other hand one can consider the Abelian projection as a gauge-dependent 
tool to associate the confining gluon configurations with the Abelian monopoles. This tool may 
work well in one gauge and may not work in another gauge. 

There are conflicting reports of the projection independence of the dual superconductor 
mechanism of the color confinement 1 . The Abelian and Monopole dominance were observed in 
more than one gauge 0^11]. The London penetration length measured in the MA projection 
is the same as the one obtained without gauge fixing [IT] . The monopole condensation studied 
with the help of a monopole creation operator was observed not only in the MA gauge of SU{2) 
gauge theory |H] but also in other gauges [12] . 

On the other hand there are indications that the monopole dynamics is affected by the choice 
of the Abelian projection. In the MA gauge the monopole trajectories percolate only in the 
confinement phase contrary to the case without any gauge fixing, in which the monopoles are 
percolating in any phase The chromoelectric string in different Abelian projections looks 
differently: the correlation length (the inverse monopole mass) extracted from the string profile 
in the AP gauge is consistent with zero contrary to the MA gauge |JJ]. The chiral condensate 
is dominated by the contributions of the Abelian monopoles in the MA gauge ^3 E] contrary 
to F12 [IH] and AP jTH] gauges. 

One can show analytically that in the AP gauge the dual superconductor mechanism can 
not be realized [9]. The reason is very simple: the monopoles in the continuum limit in this 
gauge are static and they can not contribute to the potential between static heavy quarks. On 
the other hand it was shown in Ref. [5] that the absence of the dual superconductor description 
in AP gauge does not contradict to the statement of Ref. ^2] that the monopole condensation 
is realized in any gauge. 

In this paper we are studying the effective potential for the monopole field using the Frohlich- 
Marchetti monopole creation operator. We evaluate the potential in AP and Fi 2 gauges 
and compare the results with the monopole potential obtained in the MA gauge [7j . The value 
of the monopole condensate is defined by the minimum of the effective potential. We describe 
the monopole creation operator in Section [2] and present numerical results in Section [HI Our 
conclusions are summarized in the last Section. 

A brief review of the current literature on this subject can be found in Ref. [§]• 
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2 Abelian monopole creation operator in SU(2) model 



We study the 577(2) gauge theory with the standard Wilson action, S[U] = — 1/2 J^p^r Up, 
where the sum goes over the plaquettes P and Up is the SU{2) plaquette variable composed of 
the link fields Ui, Up = UiU2U\u\. The link field is parameterized in the standard way: 




sm< 



< < vr/2 and -7r < 6, x < fl". 

In Abelian projection the residual gauge transformation matrices have the diagonal form 
fl Ahel (uj) = diag{e iuj ,e~ iuJ ), where to is an arbitrary scalar function. Under these transformations 
the diagonal field 6 transforms as an Abelian gauge field, 6 Xfl — > 6 xtl -\-uj x — u x+ jj,, the off-diagonal 
field x changes as a double charged matter field, Xxn — ► Xxp, + oj x + uj x+ fi, the field remains 
intact. The SU{2) plaquette action contains [TH] various interactions between these fields as 
well as the action for the Abelian gauge field 6: 

5[?7] = -£/3p[0]cos0p + ... . (2) 
p 

Here Op = Q\ + 6> 2 — #3 — 6> 4 is a lattice analogue of the Abelian field strength tensor and Pp[4>] 
is an effective coupling constant dependent of the fields <fi, Ref. [THj . 

Following Ref. j7j we apply the monopole creation operator of Frohlich-Marchetti to 
the Abelian part of the non- Abelian action (J2J). Effectively, this operator shifts the Abelian 
plaquette variable 6 P as follows 2 : 



^mon(x) = exp 




-cosflp + cos(#p + W P {x)) \, (3) 



where Wp = 2tt5A^ 1 (D x — u x ), uo x is a Dirac string attached to the monopole and the Dirac 
cloud D x satisfies the equation 5*D X = *5 X . 

We have used the differential form notations on the lattice described in detail in the second 
paper in Ref. [3]. S (d) is the backward (forward) derivative on the lattice which decreases 
(increases) by one the rank of the form on which it acts. The rank of the form is determined by 
a dimensionality of the lattice cell on which this form is defined. For example, a scalar function 
is a 0-form, the vector function is a 1-form etc. If A is a lattice vector, then 5 A is a scalar 
(a lattice analogue of the divergence, d^A^) while dA is an antisymmetric tensor (a lattice 
analogue of the field strength, d^A^). 6 and d operators are nilponent, i.e., 5 2 = = d 2 . 
The lattice Laplacian is A = 5d + dS, and A -1 is the inverse Laplacian. The lattice Kronecker 
symbol is denoted as 5 X : it is a scalar function which is equal to unity at the site x and zero 
otherwise. The *-operator relates the forms on the dual and original lattice. For example, if 
on the four-dimensional lattice B is a scalar function (0-form) on the original lattice, then *B 
is a 4-form on the dual lattice. 



2 Notc that in this paper we are using the "old" definition J7| of the monopole creation operator. The "new" 
definition J2] takes into account charged matter fields but it is very involved from the point of view of numerical 
calculations. Moreover, results of Ref. clearly show that there is no qualitative difference between the old 
and the new definitions. 
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The operator (j3J) is clearly gauge invariant with respect to the £7(1) gauge transformations. 
One can also perform a formal duality transformation with respect to the quantum average 
of the operator © and show that in the dual model - which has form of the Abelian Higgs 
model - this operator is invariant under the (dual) gauge transformations [7j. Moreover, one 
can represent the partition function as a sum over closed monopole trajectories. In this repre- 
sentation the quantum average of the monopole creation operator $ mon (x), Eq. (JHJ), is given by 
a sum over all closed monopole trajectories plus one open trajectory which begins at the point 
x, Refs. [17j. Thus, this operator creates a monopole at the point x. 

To get the monopole condensate we have to study the effective constraint potential for the 
monopole creation operator $ m0 n, 

V effc (<S>) = " ln«*(* " ^ E *man(x)))) . (4) 

' x 

This potential selects the zero-momentum component of the creation operator. The minimum 
of this potential corresponds to the monopole condensate. However, a numerical calculation of 
the potential V e ff c (Q) is time consuming, and in this paper we present results for the probability 
distribution 

V($) = -ln({6($-$ mon (x)))), (5) 

which has a meaning very similar to (|1J). 

We perform our study in the Abelian Polyakov gauge and in the Abelian field strength gauge 
which are defined as we already discussed by the diagonalization of the (untraced) Polyakov 
loop, -Pb[^1> an d of the £^12 plaquette, 

Px[U] = U xA U x+iA . . . U x _ lA , U Xjl2 = U Xjl U x+u Ul +%1 Ul 2 , (6) 

with respect to gauge transformations, U X4l — > f^L^fiJ. 

We compare the potential in these gauges with the monopole potential obtained in the MA 
gauge in Ref. [7j. The MA gauge is defined by the maximization of the lattice functional 

maxiW^l , Rma[U] = £ Tr (a 3 U (s, ^)a 3 U\s, //)) . (7) 

s,ji 

(03 is the Pauli matrix). In the continuum limit a local condition of the maximization can be 
written in the form of the differential equation, (<9 M + igAfy^A 1 ^ — iA 2 ^) = 0. 



3 Numerical results 

We simulate the 377(2) gauge fields on the lattices L 3 S x 4, L s = 12, 14, 16, 24 with C-periodic 
boundary conditions in space directions [2T]. The C-periodicity corresponds to the anti- 
periodicity for the Abelian gauge fields which is required by the Gauss law 3 [7j. In the case of 
SU(2) gauge group the C-periodic boundary conditions mean that on the space boundary we 
have U x ^ — > Q + U X>IX Q, Q = ia 2 - 

3 One can not introduce the creation operator of the charged particle in a finite volume with periodic boundary 
conditions. 
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We fix AP and Fu gauges as it was described in the previous Section. The results for the 
MA gauge (quoted below) are taken from Ref. |7j. To get the effective potential in the AP 
and F 12 gauges we use 400 independent configurations of SU (2) gauge fields for each value 
of the gauge coupling (3 at a given lattice volume. On each configuration the distribution of 
the monopole creation operator is evaluated in 20 points. The logarithm of the distribution 
provides us with the effective potential 

To evaluate the numerical errors of the potential we use the bootstrap method. For each 
value of (3 and the lattice volume we get a distribution of the values of the monopole creation 
operator ("initial ensemble"). Then we use the initial ensemble to construct a number (typically 
N ens = 500) of additional ensembles of the operator values by randomly choosing the entries 
from the initial ensemble. Each value from the initial ensemble may enter the additional 
ensemble more than one time. The number of entries in each of the additional ensembles is 
the same as in the initial one. Then for each ensemble we construct a histogram, the (minus) 
logarithm of which has a meaning of the monopole potential, V($), according to Eq. (jSJ). Thus, 
for each value of the monopole field, $, we get N ens values of the potential, V, which form a 
Gaussian distribution. The central value of this distribution gives us the value of the potential 
at given lattice volume L 2 S x L t , (3 and $, V = V(Q), while the width of the distribution is the 
statistical error. 

The examples of the effective monopole potentials for 16 3 x 4 lattice are shown on Figs.^a- 
d). We depict the positive part of the potentials ($ > 0) at various values of the gauge 
coupling (3 in the MA, AP and F\2 gauges. The minimum of the effective potential corresponds 
to the value of the monopole condensate. The critical gauge coupling corresponding to the 
temperature phase transition at our lattices is [3 C ~ 2.3. Thus Figs. ma,b,c) correspond to the 
confinement phase while Fig. ^d) corresponds to the deconfinement phase. 

First of all we note that for all considered values of (3 (i) the minima of the potentials in 
the AP and F 12 gauges coincide with each other within numerical errors; (ii) the potential in 
the MA gauge is different from AP and F\i potentials. According to Fig. E^a) in the strong 
coupling limit (j3 = 0.1) the minima of the monopole potential in all three gauges are located at 
the same point, $ m i n ~ 1. As we increase (3 the difference in the monopole condensates in MA 
gauge and in AP and F i2 gauges appears, see Fig. Hfb,c). Moreover, in the deep deconfinement 
phase, Fig. E^d), the monopole condensate vanishes in MA gauge while in AP and F 12 gauges 
the condensate is non-zero. 

Since the phase transition in the SU(2) gauge theory is of the second order, the finite volume 
effects may be essential for the determination of the monopole condensate. To get rid of finite 
volume corrections we measure the condensate on the lattices with various spatial extensions 
(L s = 12, 14, 16, 24) and extrapolate the value of the condensate to the thermodynamic limit, 
L s — > oo: 

$ c = V? + j. (8) 

The examples of the fits for the AP and the F 12 gauges are shown in Figs. |2fa,b). The values 
of x 2 /d.o.f. are in the range 0.2 ~ 1. 

The monopole condensates in the thermodynamic limit (L s — > oo) for all three Abelian 
projections are shown in Fig. |3] as functions of (3. One can clearly see that the monopole 
condensate in the MA projection vanishes at a certain critical (3 = (3 C which is very close to 
the phase transition point, (3 C ~ 2.3. Contrary to the MA gauge the monopole condensates 



5 



: 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r-a c — i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 



(c) (d) 

Figure 1: The potential on the monopole field in the MA, AP and F 12 gauges at various values 
of the gauge coupling (3. The data for the MA gauge are taken from Ref. [Zj. 

obtained in the AP and the F\2 gauges do not vanish at (5 — (3 C . This result is in contradiction 
with observations of Ref. [T2] . 

The dependence of the monopole condensate on (3 can be fitted by the following function: 

&?{{3) = l-{(3/f3 c )\ (9) 

with x 2 /d.o.f. ~ 0.3. It occurs that 7 = 1.2(5) and /3 C = 2.31(3). The value of (3 C coincides 
within error bars with the known critical value [22| on x 4 lattice. 
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Figure 2: The extrapolation (JHJ) of the monopole condensate to the thermodynamic limit in the 
AP and the gauges for various values of the gauge coupling (3. The spatial extensions of 
the lattices are L s = 12, 14, 16, 24. 
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Figure 3: The monopole condensate in the thermodynamic limit in the MA, the AP and the 
F12 gauges. The dash-dotted line is the fit of the monopole condensate in the MA gauge by 
Eq. ((0]). The critical value of the gauge coupling (along with the numerical error) is denoted 
by the vertical dashed line. 
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4 Conclusion 



Summarizing, we have presented an evidence that the monopole condensate in different Abelian 
projections coincide with each other only in the (unphysical) strong coupling region. Generally, 
the condensate depends on the choice of the Abelian projection. We have considered three 
Abelian projections and only in MA projection the vacuum behaves as the dual superconductor. 
Our results are in contradiction with conclusions of Ref. jT^j where condensate was found to 
be projection-independent. 
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